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Abstract. The positive and not completely positive maps of density matrices, 
which are contractive maps, are discussed as elements of a semigroup. A new kind 
of positive map (the purification map), which is nonlinear map, is introduced. 
The density matrices are considered as vectors, linear maps among matrices are 
represented by superoperators given in the form of higher dimensional matrices. 
Probability representation of spin states (spin tomography) is reviewed and U{N)- 
tomogram of spin states is presented. Properties of the tomograms as probability 
distribution functions are studied. Notion of tomographic purity of spin states is 
introduced. Entanglement and separability of density matrices are expressed in 
terms of properties of the tomographic joint probability distributions of random 
spin projections which depend also on unitary group parameters. A new positivity 
criterion for hermitian matrices is formulated. An entanglement criterion is given 
in terms of a function depending on unitary group parameters and semigroup of 
positive map parameters. The function is constructed as sum of moduli of U{N)- 
tomographic symbols of the hermitian matrix obtained after action on the density 
matrix of composite system by a positive but not completely positive map of the 
subsystem density matrix. Some two-qubit and two-qutritt states are considered as 
examples of entangled states. The connection with the star-product quantisation is 
discussed. The structure of the set of density matrices and their relation to unitary 
group and Lie algebra of the unitary group are studied. Nonlinear quantum evolution 
of state vector obtained by means of applying purification rule of density matrices 
evolving via dynamical maps is considered. Some connection of positive maps and 
entanglement with random matrices is discussed and used. 

Keywords: unitary group, entanglement, adjoint representation, tomogram, oper- 
ator symbol, random matrix. 



1. Introduction 

The states in quantum mechanics are associated with vectors in Hilbert 
space [1] (it is better to say with rays) in the case of pure states. 
For mixed state, one associates the state with density matrix [2, 3]. 
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In classical mechanics (statistical mechanics), the states are associ- 
ated with joint probability distributions in phase space. There is an 
essential difference in the concept of states in classical and quantum 
mechanics. This difference is clearly pointed out by the phenomenon of 
entanglement. The notion of entanglement [4] is related to the quantum 
superposition rule of the states of subsystems for a given multipartite 
system. For pure states, the notion of entanglement and separability 
can be given as follows. 

If the wave function [5] of a state of a bipartite system is represented 
as the product of two wave functions depending on coordinates of the 
subsystems, the state is simply separable; otherwise, the state is sim- 
ply entangled. An intrinsic approach to the entanglement measure was 
suggested in [6]. The measure was introduced as the distance between 
the system density matrix and the tensor product of the associated 
states. For the subsystems, the association being realized via partial 
traces. There are several other different characteristics and measures 
of entanglement considered by several authors [7-13]. For example, 
there are measures related to entropy (see, [14-24]). Also linear entropy 
of entanglement was used in [25-27], "concurrences" in [28, 29] and 
"covariance entanglement measure" in [30] . Each of the entanglement 
measures describes some degree of correlation between the subsystems' 
properties. 

The notion of entanglement is not an absolute notion for a given 

system but depends on the decomposition into subsystems. The same 
quantum state can be considered as entangled, if one kind of division 
of the system into subsystems is given, or as completely disentangled, 
if another decomposition of the system into subsystems is considered. 

For instance, the state of two continuous quadratures can be entan- 
gled in Cartesian coordinates and disentangled in polar coordinates. 
Coordinates are considered as measurable observables labelling the 
subsystems of the given system. The choice of different subsystems 
mathematically implies the existence of two different sets of the sub- 
systems' characteristics (we focus on bipartite case). We may consider 
the Hilbert space of states H{1,2) or H{1',2'). The Hilbert space for 
the total system is, of course, the same but the index (1,2) means 
that there are two sets of operators Pi and P2, which select subsystem 
states 1 and 2. The index (1', 2') means that there are two other sets 
of operators P{ and P2, which select subsystem states 1' and 2'. The 
operators Pi^2 and P[, 2/ have specific properties. They are represented 
as tensor products of operators acting in the space of states of the 
subsystem 1 (or 2) and unit operators acting in the subsystem 2 (or 
1). In other words, we consider the space H, which can be treated as 
the tensor product of spaces H(l) and H{2) or H(l') and H{2'). In 
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the subsystems 1 and 2, there are basis vectors | ni) and | 1112), and in 
the subsystems 1' and 2' there arc basis vectors | n'l) and | m'-j). The 
vectors | ni) | 7712) and | n'^^) | 77X2) form the sets of basis vectors in 
the composite Hilbert space, respectively. These two sets are related 
by means of unitary transformation. An example of such a composite 
system is a bipartite spin system. 

If one has spin-ji [the space H{1)] and spin-_72 [the space H{2)] 
systems, the combined system can be treated as having basis 
I hmi) I J2?n2). 

Another basis in the compositc-systcm-state space can be considered 
in the form | jm) , where j is one of the numbers | ji — j2 1 , | ji — j2 1 + 1 , • • • , 
Ji + J2 and m = m\ + m2- The basis | jm) is related to the basis 
I ji^i) I i2?7i2) by means of the unitary transform given by Clcbsch- 
Gordon coefficients C {jimij2m2\jm). From the viewpoint of the given 
definition, the states | jm) are entangled states in the original basis. 
Another example is the separation of the hydrogen atom in terms of 
parabolic coordinates used while discussing the Stark effect. 

The spin states can be described by means of the tomographic 
map [31-33]. For bipartite spin systems, the states were described 
by the tomographic probabilities in [34, 35]. Some properties of the 
tomographic spin description were studied in [36]. In the tomographic 
approach, the problems of the quantum state entanglement can be 
cast into the form of some relations among the probability distribution 
functions. On the other hand, to have a clear picture of entanglement, 
one needs a mathematical formulation of the properties of the density 
matrix of the composite system, a description of the linear space of 
the composite system states. Since a density matrix is hermitian, the 
space of states may be embedded as a subset of the Lie algebra of 
the unitary group, carrying the adjoint representation of U{p?), where 
v? = (2j + 1)^ is the dimension of the spin states of two spinning 
particles. Thus one may try to characterize the entanglement phenom- 
ena by using various structures present in the space of the adjoint 
representation of the U{n'^) group. 

The aim of this paper is to give a review of different aspects of 
density matrices and positive maps and connect entanglement prob- 
lems with the properties of tomographic probability distributions and 
discuss the properties of the convex set of positive states for composite 
system by taking into account the subsystem structures. We used [6] the 
Hilbert-Schmidt distance to calculate the measure of entanglement as 
the distance between a given state and the tensor product of the partial 
traces of the density matrix of the given state. In [37] another measure 
of entanglement as a characteristic of subsystem correlations was intro- 
duced. This measure is determined via the covariance matrix of some 
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observables. A review of different approaches to the entanglement no- 
tion and entanglement measures is given in [38] where the approach to 
describe entanglement and separability of composite systems is based, 
e.g., on entropy methods. 

Due to a variety of approaches to the entanglement problem, one 
needs to understand better what in reality the word "entanglement" 
describes. Is it a synonym of the word "correlation" between two sub- 
systems or does it have to capture some specific correlations attributed 
completely and only to the quantum domain? 

The paper is organized as follows. 

In section 2 we discuss division of composite systems onto subsys- 
tems and relation of the density matrix to adjoint representation of 
unitary group in generic terms of vector representation of matrices; 
wc study also completely positive maps of density matrices. In sec- 
tion 3 we consider a vector representation of probability distribution 
functions and notion of distance between the probability distributions 
and density matrices. In section 4 we present definition of separable 
quantum state of a composite system and criterion of separability. In 
section 5 the entanglement is considered in terms of operator symbols. 
Particular tomographic probability representation of quantum states 
and tomographic symbols are reviewed in section 6. Symbols of mul- 
tipartite states are studied in section 7. In section 8 spin tomography 
is reviewed. An example of qubit state is done in section 9. The uni- 
tary spin tomogram is introduced in section 10 while in section 11 
dynamical map and corresponding quantum evolution equations are 
discussed as well as examples of concrete positive maps. Conclusions 
and perspectives are presented in section 12. 



2. Composite system 

In this section, we review the meaning and notion of composite system 
in terms of additional structures on the linear space of state for the 
composite system. 

2.1. Difference of states and observables 

In quantum mechanics, there are two basic aspects, which are associ- 
ated with linear operators acting in a Hilbert space. The first one is 
related to the concept of quantum state and the second one, to the 
concept of observable. These two concepts of state and observable are 
paired via a map with values in probability measures on the real line. 
Often states are described by Hermitian nonnegative, trace-class, ma- 
trices. The observables are described by Hermitian operators. Though 
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both states and observables are identified with Hermitian objects, there 
is an essential difference between the corresponding objects. The ob- 
servables have an additional product structure. Thus we may consider 
the product of two linear operators corresponding to observables. 

For the states, the notion of product is redundant. The product of 
two states is not a state. For states, one keeps only the linear structure 
of vector space. For finite n-dimensional system, the Hermitian states 
and the Hermitian observables may be mapped into the Lie algebra 
of the unitary group U{n). But the states correspond to nonnegative 
Hermitian operators. Observables can be associated with both types 
of operators, including nonnegative and nonpositive ones. The space 
of states is therefore a convex-linear space which, in principle, is not 
equipped with a product structure. Due to this, transformations in the 
linear space of states need not preserve any product structure. In the set 
of observables, one has to be concerned with what is happening with the 
product of operators when some transformations are performed. State 
vectors can be transformed into other state vectors. Density operators 
also can be transformed. We will consider linear transformations of the 
density operators. The density operator has nonnegative eigenvalues. In 
any representation, diagonal elements of density matrix have physical 
meaning of probability distribution function. 

Density operator can be decomposed as a sum of eigenprojectors 
with coefficients equal to its eigenvalues. Each one of the projectors 
defines a pure state. There exists a basis in which every eigenprojector 
of rank one is represented by a diagonal matrix of rank one with only 
one matrix element equal to one and all other matrix elements equal 
to zero. Other density matrices with similar properties belong to the 
orbit of the unitary group on the starting eigenprojector. Depending on 
the number of distinct nonzero values determines the class of the orbit. 
Since density matrices of higher rank belong to an appropriate orbit 
of a convex sum of the different diagonal eigenprojectors (in special 
basis), we may say that generic density matrices belong to the orbits of 
the unitary group acting on the diagonal density matrices which belong 
to the Cartan subalgebra of the Lie algebra of the unitary group. Any 
convex sum of density matrices can be treated as a mean value of a 
random density matrix. The positive coefficients of the convex sum can 
be interpreted as a probability distribution function which makes the 
averaging providing the final value of the convex sum. The set of density 
matrices may be identified with the union of the orbits of the unitary 
group acting on diagonal density matrices considered as elements of 
the Cartan subalgebra. 
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2.2. Matrices as vectors, density operators and 
superoperators 

When matrices represent states it may be convenient to identify them 
with vectors. In this density matrix can be considered as a vector 

with additional properties of its components. If the identifications are 
done elegantly, we can see the real Hilbert space of density matrices in 
terms of vectors with real components. In this case, linear transforms 
of the matrix can be interpreted as matrices called superoperators. It 
means that density matrices-vectors undergoing real linear transfor- 
mations are acted on by the matrices representing the action of the 
superoperators of the linear map. This construction can be continued. 
Thus we can get a chain of vector spaces of higher and higher dimen- 
sions. Let us first introduce some extra constructions of the map of 
a matrix onto a vector. Given a rectangular matrix M with elements 
Mid, where i = 1,2, ... ,n and d = 1,2, ... ,m, one can consider the 
matrix as a vector A4 with N = nm components constructed by the 
following rule: 

Ml = Mil, M2 = M12, Mm = Mim, 

(1) 

Mm+l = M21, . . . Mn = Mnm- 

Thus we construct the map M ^ M =tj^j^M. 

We have introduced the linear operator t^^j^ which maps the matrix 

M onto a vector M. Now we introduce the inverse operator Pj^j^,f which 
maps a given column vector in the space with dimension = mn onto a 
rectangular matrix. This means that given a vector M = M\, . . . ,Mn, 
wc relabel its components by introducing two indices i = 1, . . . ,n and 
d = l,...,m. The relabeling is accomplished according to (1). Then 
we collect the relabeled components into a matrix table. Eventually we 
get the map 

PmmM = M. (2) 
The composition of these two maps 

iuMPMM^ =1 • ^ (3) 

acts as the unit operator in the linear space of vectors. 

Given a nxn matrix the map considered can also be applied. The 
matrix can be treated as an n^-dimensional vector and, vice versa, the 
vector of dimension may be mapped by this procedure onto the nxn 
matrix. 
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Let us consider a linear operator acting on the vector M. and related 
to a linear transform of the matrix M. First, we study the correspon- 
dence of the linear transform of the form 

M ^ gM = Ml (4) 

to the transform of the vector 

M^M^g = C^gM. (5) 

One can show that the x matrix £^ is determined by the tensor 
product of the nxn matrix g and nxn unit matrix, i.e., 

4 = 5(^1. (6) 
Analogously, the linear transform of the matrix M of the form 

M ^Mg = Ml (7) 

induces the linear transform of the vector M of the form 

where the x matrix Cg reads 

Cl = l^g^'. (9) 
Similarity transformation of the matrix M of the form 

M ^ gMg-^ (10) 
induces the corresponding linear transform of the vector Ai of the form 

M^Ms=ClM, (11) 
where the n? x matrix Cg reads 

>C^ = 5®(5-')*'- (12) 

Starting with vectors, one may ask how to identify on them a product 
structure which would make p into an algebra homomorphism. An 
associative algebraic structure on the vector space may be defined by 
imitating the procedure one uses to define star-products on the space 
of functions on phase space. One can define the associative product of 
two A'^- vectors Mi and M2 using the rule 

M = Mii^M2, (13) 
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where 

N 

Mk=Yl Kt{Mi)i{M2)s. (14) 

l,s=l 

If one applies a linear transform to the vectors Mi, M.2, M. oi the form 

Mi^M[ = CMi, M2^M2 = CM2, M ^ M' = CM, 

and requires the invariance of the star-product kernel, one finds 

M'iicM'2=M', if >C = G(8)G-^*^ GeGL{n). 

The kernel K^^ (structure constants) which determines the associative 
star-product satisfies a quadratic equation. Thus if one wants to make 
the correspondence of the vector star-product to the standard matrix 
product (row by column), the matrix M must be constructed appro- 
priately. For example, if the vector star-product is commutative, the 
matrix M corresponding to the A'^- vector M can be chosen as a diagonal 
NxN matrix. This consideration shows that the map of matrices on the 
vectors provides the star-product of the vectors (defining the structure 
constants or the kernel of the star-product) and, conversely, if one starts 
with vectors and uses matrices with the standard multiplication rule, 
it will be the map to be determined by the structure constants (or by 
the kernel of the vector star-product). 

The constructed space of matrices associated with vectors enables 
one to enlarge the dimensionality of the group acting in the linear space 
of matrices in comparison with the standard one, i.e., we may relax the 
requirement of invariance of the product structure. In general, given 
a nxn matrix M the left action, the right action, and the similarity 
transformation of the matrix are related to the complex group GL{n). 
On the other hand, the linear transformations in the linear space of v?- 
vectors M obtained by using the introduced map are determined by 
the matrices belonging to the group GL{'n?). There are transformations 
on the vectors which cannot be simply represented on matrices. If 
M — > $(M) is a linear homogeneous function of the matrix M, we 
may represent it by 

^ab = Baa',bb'Ma'b'- 

Under rather clear conditions, Baa',bb' can be expressed in terms of its 
nonnormalized left and right eigenvectors: 

V 

being an index for eigenvalues, which corresponds to 
$(M) = xMy^ = ^ x{y)My\v). 
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There are possible linear transforms on the matrices and correspond- 
ing hnear transforms on the induced vector space which do not give rise 
to a group structure but possess only the structure of algebra. One can 
describe the map of nxn matrices M (source space) onto vectors Ai 
(target space) using specific basis in the space of the matrices. The basis 
is given by the matrices Ejk (j, k = 1,2, ... ,n) with all matrix elements 
equal to zero except the element in the j'th row and kth column which 
is equal to unity. One has the obvious property 

Mjk = Tr (MEjk) . (15) 

In our procedure, the basis matrix Ejk is mapped onto the basis column- 
vector Sjk, which has all components equal to zero except the unity 
component related to the position in the matrix determined by the 
numbers j and k. Then one has 

n 

M=Y1 Ti (MEjk) £jk. (16) 
j,k=i 

For example, for similarity transformation of the finite matrix M, one 
has 

N 

•^9 = E Tr {gMg-^Ejk) Sjk- (17) 

j,k=i 

Now we will define the notion of 'composite' vector which corre- 
sponds to dividing a quantum system into subsystems. 
We will use the following terminology. 

In general, the given linear space of dimensionality = mn has 
a structure of a bipartite system, if the space is equipped with the 
operator Pj^j^ and the matrix M (obtained by means of the map) has 
matrix elements in factorizable form 

Mid XiVd. (18) 

This M = x®y corresponds to the special case of nonentangled states. 
Otherwise, one needs 

M = ^x{v) (g) y{u). 

In fact, to consider in detail the entanglement phenomenon, in the 
bipartite system of spin- 1/2, one has to introduce a hierarchy of three 
linear spaces. The first space of pure spin states is the two-dimensional 
linear space of complex vectors 



.Xl 
X2 



(19) 
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In this space, the scalar product is defined as follows: 

{x\y} = xlvi + x*2y2. 



(20) 



So it is a two-dimensional Hilbcrt space. Wc do not equip this space 
with a vector star-product structure. In the primary linear space, one 
introduces linear operators M which are described by 2x2 matrices 
M. Due to the map discussed in the previous section, the matrices 
are represented by 4-vectors M. belonging to the second complex 4- 
dimensional space. The star-product of the vectors M. determined by 
the kernel JCf^ is defined in such a manner in order to correspond to 
the standard rule of multiplication of the matrices. 

In addition to the star-product structure, we introduce the scalar 
product of the vectors Mi and M2, in view of the definition 



{Ml I M2) = Ti{MlM2) 



(21) 



which is the trace formula for the scalar product of matrices. 

This means introducing the real metric g'^^ in the standard notation 
for scalar product 



{Ml \M2)= E iMi)l9''^iM2)/3, 

a,f3=l 



(22) 



where the matrix of the form 



/ 1 \ 

10 

10 

VO 1/ 



(23) 



The scalar product is invariant under the action of the group of non- 
singular 4x4 matrices £, which satisfy the condition 



9 



(24) 



The product of matrices £ satisfies the same condition since g'^ = 1. 

Thus, the space of operators M in the primary two-dimensional 
space of spin states is mapped onto the linear space which is equipped 
with a scalar product (metric Hilbert space structure) and an associa- 
tive star-product (kernel satisfying the quadratic associativity equa- 
tion). In the linear space of the 4-vectors M, we introduce linear 
operators (superoperators), which can be associated with the algebra 
of 4x4 complex matrices. 

Let us now focus on density matrices. This means that our matrix M 
is considered as a density matrix p which describes a quantum state. 
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We consider here the action of the unitary transformation U{n) of 
the density matrices and corresponding transformations on the vector 
space. If one has the structure of a bipartite system, we also consider 
the action of local gauge transformation both in the "source space" of 
density matrices and in the "target space" of the corresponding vectors. 
The nxn density matrix p has matrix elements 

Pik = ph, Trp = l, (VHV')>0. (25) 

Since the density matrix is hermitian, it can always be identified as an 
element of the convex subset of the linear space associated with the Lie 
algebra of U{n) group, on which the group U{n) acts with the adjoint 
representation 

p^pu = UpUl (26) 

The system is said to be bipartite if the space of representation is 
equipped with an additional structure. This means that for 

= rii ■ 

where, for simplicity, rii = n2 = n, one can make first the map of 
nxn matrix p onto n^-dimensional vector p according to the previous 
procedure, i.e., one equips the space by an operator t^-p. Given this 
vector one makes a relabeling of the vector p components according to 
the rule 

P^Pid,ke^ i,k = 1,2, ... ,ni, ci,e = 1,2, ... ,77,2, (27) 
i.e., obtaining again the quadratic matrix 

Pq = Pp,pP- (28) 

The unitary transform (26) of the density matrix induces a linear 
transform of the vector p of the form 

p^pu = {U®U*)p. (29) 

There exist linear transforms (called positive maps) of the density ma- 
trix, which preserve its trace, hermicity, and positivity. In some cases, 
they have the following form introduced in [39] 

Po^ Pu = Lupo = ^PkUkPoUl, 'YjPk = 1' (30) 

k k 

where are unitary matrices and are positive numbers. 

If the initial density matrix is diagonal, i.e., it belongs to the Cartan 
subalgebra of the Lie algebra of the unitary group, the diagonal ele- 
ments of the obtained matrix give a "smoother" probability distribution 
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than the initial one. A generic transformation preserving previously 
stated properties may be given in the form (see [39, 40]) 

Po^Pv = Lvpo = E VkPoV^, Yl ^k^k = 1- (31) 

k k 

For example, if Vk {k = 1, 2, . . . ,N) are taken as square roots of or- 
thogonal projectors onto complete set of N state, the map provides 
the map of the density matrix po onto diagonal density matrix pod 
which has the same diagonal elements as po has. In this case, the 
matrices Vk have the only nonzero matrix element which is equal to 
one. Such a map may be called "decoherence map" because it removes 
all nondiagonal elements of the density matrix po killing any phase 
relations. In quantum information terminology, one uses also the name 
"phase damping channel." More general map may be given if one takes 
Vk as N generic diagonal density matrices, in which eigenvalues are 
obtained by N circular permutations from the initial one. Due to this 
map, one has a new matrix with the same diagonal matrix elements 
but with changed nondiagonal elements. The purity of this matrix is 
smaller then the purity of the initial one. This means that the map 
is contractive. All matrices with the same diagonal elements up to 
permutations belong to a given orbit of the unitary group. 

For a large number of terms with randomly chosen matrices Vk in 
the sum in (31), the above map gives the most stochastic density matrix 

Po^ Ps = Lipo = (n)"h. 

Its four-dimensional matrix Li for the qubit case has four matrix el- 
ements different from zero. These matrix elements are equal to one. 
They have the labels Ln, L14, L41, L44. The map with two nonzero 
matrix elements L41 = L44 = provides pure-state density matrix from 
any po- The transform (30) is the partial case of the transform (31). We 
discuss the transforms separately since they are used in the literature 
in the presented form. 

One can see that the constructed map of density matrices onto 
vectors provides the corresponding transforms of the vectors, i.e., 

Po^ Pu = Ypk{Uk®Uk)pQ (32) 

k 

and 

P^^Pv = Y.i^k®^kW (33) 

k 

It is obvious that the set of linear transforms of vectors, which preserve 
their properties of being image of density matrices, is essentially larger 
than the standard unitary transform of the density matrices. 
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Formulae (32) and (33) mean that the positive map superoperators 
acting on the density matrix in the vector representation are described 
by X matrices 

k 



and 

k 

respectively. 

The positive map is called "noncompletely positive" if 



(35) 



k 



EnV.-E 



vlv. 



This map is related to a possible "nonphysical" evolution of a subsys- 
tem. 

Formula (34) can be considered in the context of random matrix rep- 
resentation. In fact, the matrix Cu can be interpreted as the weighted 
mean value of the random matrix Uk U^. The dependence of matrix 
elements and positive numbers pk on index k means that we have a 
probability distribution function pk and averaging of the random matrix 
Uk (g) by means of the distribution function. So the matrix jCu reads 



Cu = {U(^U*). 



(36) 



Let us consider an example of a 2x2 unitary matrix. We can consider 
a matrix of the SU{2) group of the form 



u 



a j3 
-P* a* 



(37) 



The 4x4 matrix Ctt takes the form 



f £ m 
—n s 
-n* —q* 

-m - 



m 

-Q 

s* 



n 



n 



(38) 



m 



The matrix elements of the matrix Cu are the means 

m = (a/?*), 

e = {aa*), 

n = {af5), 

s = (a^), 

9 = iP'). 



(39) 
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The moduli of these matrix elements are smaller than unity. 
The determinant of the matrix Cu reads 



deiCu = (1 - 2i){\q 



+ 4Re 



q m n + runs 



If one represents the matrix Cu in block form 



A B 
C D 



then 



and 



A 



I m 
—n s 



D = (72^V2, 



B 



m* 1 - 
—q n 



C = -a2B*a2, 



(40) 



(41) 



(42) 



(43) 



where o"2 is the Pauli matrix. 

One can check that the product of two different matrices Cu can 
be cast in the same form. This means that the matrices Cu form a 
9-parameter compact semigroup. It means that the product of two 
matrices from the set (semigroup) belongs to the same set. It means 
that composition is inner like the one for groups. There is a unity 
element in the semigroup, however, there exist elements which have no 
inverse. In our case, these elements arc described, e.g., by the matrices 
with zero determinant. Also the elements, which are matrices with 
nonzero determinants, have no inverse elements in this set, since the 
map corresponding to the inverse of these matrices is not positive one. 
For example, in the case £ = 1/2 and m = 0, one has the matrices 



A 



1/2 
— n s 



B 



1/2 
—q n 



(44) 



The determinant of the matrix Cjj in this case is equal to zero. All the 
matrices Cjj have the eigenvector 



/l/2\ 



Vl/2/ 



I.e., 



^uPo = Po- 

This eigenvector corresponds to the density matrix 



Pi 



1/2 
1/2 



(45) 



(46) 



(47) 
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which is obviously invariant of the positive map. 

For random matrix, one has correlations of the random matrix ele- 
ments, e.g., {aa*) ^ (a) (a*). 



The matrix Cr, 



/I 




Vo 





1 

1 







1/ 



(48) 



maps the vector 



/Pii\ 

Pl2 
P21 
\ P22 J 



(49) 



onto the vector 



Pt 



(50) 



P21 
Pl2 
\ P22 I 

This means that the positive map (48) connects the positive density 
matrix with its transpose (or complex conjugate). This map can be 
presented as the connection of the matrix p with its transpose of the 
form 



„tr 



\{p-^ 0-lPCri - 0-2p(72 + (73/9(73 j . 



There is no unitary transform connecting these matrices. 

There is noncompletely positive map in the A'^-dimensional case, 
which is given by the generalized formula (for some e) 

P^ Ps = -£p + -j^ Iat. 

In quantum information terminology, it is called "depolarizing chan- 
nel." 

For the qubit case, matrix form of this map reads 



L 



f 


1-6 

2 








l+£ 
2 







— e 
















— e 





V 


l+e 
2 








1-e 
2 



(51) 



Thus we constructed the matrix representation of the positive map of 
density operators of the spin- 1/2 system. This particular set of matrices 
realize the representation of the semigroup of real numbers — 1 < e < 1. 
If one considers the product £i£2 = £3, the result £3 belongs to the 
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semigroup. Only two elements 1 and —1 of the semigroup have the 
inverse. These two elements form the finite subgroup of the semigroup. 
The semigroup itself without element e = can be embedded into the 
group of real numbers with natural multiplication rule. Each matrix 
L has an inverse element in this group but all the parameters of the 
inverse elements rj live out of the segment —1,1. The group of the real 
numbers is commutative. The matrices of the nonunitary representation 
of this group commute too. It means that we have nonunitary reducible 
representation of the semigroup which is also commutative. To con- 
struct this representation, one needs to use the map of matrices on the 
vectors discussed in the previous section. Formulae (31) and (35) can 
be interpreted also in the context of the random matrix representation, 
but we use the uniform distribution for averaging in this case. So one 
has equality (35) in the form 

Cy = {V^ V*) (52) 

and the equality 

(Fty) = 1, (53) 

which provides constraints for the random matrices V used. 

Using the random matrix formalism, the positive (but not com- 
pletely positive) maps can be presented in the form 

C = {V ^jV*) - {v0v*), (V^V) - {vK) = 1. 

One can characterize the action of positive map on a density matrix p 
by the parameter 

. = 5^ = ^<1. 

As a remark we note that in [39] the positive maps (30) and (31) were 
used to describe the non-Hamiltonian evolution of quantum states for 
open systems. 

We have to point out that, in general, such evolution is not described 
by first-order-in-time differential equation. As in the previous case, if 
there are additional structures for the matrix in the form 

Pid,ke XiUdZkie, (54) 

which means associating with the initial linear space two extra lin- 
ear spaces where Xi,zi. are considered as vector components in the 
ni-dimensional linear space and ya, tg are vector components in n2- 
dimensional vector space, we see that one has bipartite structure of 
the initial space of state [bipartite structure of the space of adjoint 
representations of the group U{n)]. 
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Usually the adjoint representation of any group is defined per se 
without any reference to possible substructures. Here we introduce 
the space with extra structure. In addition to being the space of the 
adjoint representation of the group U (n), it has the structure of a bipar- 
tite system. The generalization to multipartite (A/"-partite) structure is 
straightforward. One needs only the representation of positive integer 

in the form 

N 

n^=Y[ 4. (55) 
k=l 

If one considers the more general map given by superoperator (35) 
rewritten in the form 

Cy = {V®V*), (fV)=i, 

the number of parameters determining the matrix Cy can be easily 
evaluated. For example, for n = 2, 

V - f ^ ^\ - ( ^* 

^ ~\c dj' ^ ~\c*d* 

where the matrix elements are complex numbers, the normalization 
condition provides four constraints for the real and imaginary parts of 
the matrix elements of the following matrix: 

/(|a|2) {ah*) {ha*) {bb*) \ 

{ac*) {ad*) {he*) {bd*' 

{ca*) {cb*) {da*) {db* 

V (cc*) {cd*) {dc*) {dd*) ) 



namely, 

(|a|2) + (|c|2) = 1, + = 1, {a*h) + {c*d) = 0. 

Due to the structure of the matrix £y , there are six complex parameters 

{ah*), {ac*), {ad*), {be*), {bd*), {cd*) 

or 12 real parameters. 

The geometrical picture of the positive map can be clarified if one 
considers the transform of the positive density matrix into another 
density matrix as the transform of an ellipsoid into another ellipsoid. A 
generic positive transform means a generic transform of the ellipsoid, 
which changes its orientation, values of semiaxis, and position in the 
space. But the transform does not change the ellipsoidal surface into a 
hyperboloidal or paraboloidal surface. For pure states, the positive den- 
sity matrix defines the quadratic form which is maximally degenerated. 
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In this sense, the "eUipsoid" includes all its degenerate forms corre- 
sponding to the density matrix of rank less than n (in n-dimensional 
case). The number of parameters defining the map {V (8) V*) in the 
n-dimensional case is equal to n?{n^ ~ !)• 

The linear space of Hermitian matrices is also equipped with the 
commutator structure defining the Lie algebra of the group U{n). The 
kernel that defines this structure (Lie product structure) is determined 
by the kernel that determines the star-product. 

3. Distributions as vectors 

In quantum mechanics, one needs the concept of distance between the 
quantum states. In this section, we consider the notion of distance 
between the quantum states in terms of vectors. First, let us discuss 
the notion of distance between conventional probability distributions. 
This notion is well known in the classical probability theory. 

Given the probability distribution P(k), k = 1,2, ...iV, one can 
introduce the vector P in the form of a column with components Pi = 
P{1), P2 = P(2), Pn = P{N). The vector satisfies the condition 

N 

E^fe = l- (56) 
k=l 

This set of vectors does not form a linear space but only a convex sub- 
set. Nevertheless, in this set one can introduce a distance between two 
distributions by using the one suggested by the vector space structure 

of the ambient space: 

L»2 = (p, - ^2) ' = E P^kPik + E P'^kP2k - 2 E PikP2k. (57) 

k k k 

Of course, one may use other identifications of distributions with 
vectors. 

Since all P{k) > 0, one can introduce Vk = \/P{k) as components 
of the vector V. The V can be thought of as a column with nonnegative 
components. Then the distance between the two distributions takes the 
form 

p2 = (f^ _ = 2 -2j2^/Pi{k)P2{k). (58) 

k 

The two different definitions (56) and (57) can be used as distances 
between the distributions. 

Let us discuss now the notion of distance between the quantum 
states determined by density matrices. In the density-matrix space (in 
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the set of linear space of the adjoint U{n) representation), one can 
introduce distances analogously. The first case is 

Tr{p^-P2)^ = D^ (59) 

and the second case is 

T^iVp^-Vp^f = 'D\ (60) 

In fact, the distances introduced can be written naturally as norms of 
vectors associated to density matrices 

= |pi - P2f (61) 

and 

^'=((v^i)-(v^2))'> (62) 

respectively. 

In the above expressions, wc use scalar product of vectors pi and p2 
as well as scalar products of vectors (y^^) and {^^), respectively. 

Both definitions immediately follow by identification of either matri- 
ces pi and p2 with vectors according to the map of the previous sections 

or matrices and with vectors. Since the density matrices pi 
and p2 have nonncgative eigenvalues, the matrices and y^ are 

defined without ambiguity. This means that the vectors {^/p-^) and 
iVP2^ are also defined without ambiguity. It is obvious that using this 
construction and introducing linear map of positive vectors one in- 
duces nonlinear map of density matrices. Other analogous functions, in 
addition to square root function, can be used to create other nonlinear 
positive maps. 



4. Separable systems and separability criterion 

According to the definition, the system density matrix is called sep- 
arable (for composite system) but not simply separable, if there is 
decomposition of the form 

PAB = ^Pk{pf^P^B^), E^^fc = l' l>Pfe>0. (63) 

k k 

This is Hubert's problem of biquadratcs. Is a positive biquadratic the 
positive sum of products of positive quadratics? In this formula, one 
may use also sum over two different indices. Using another labelling in 
such sum over two different indices, this sum can be always represented 
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as the sum over only one index. The formula does not demand orthogo- 
nality of the density operators p"^ and p"^ for different k. Since every 
density matrix is a convex sum of pure density matrices, one could 

demand that and p^'' be pure. This formula can be interpreted in 
the context of random matrix representation reading 

PAB = {pA® Pb), (64) 

where pA and ps are considered as random density matrices of the 
subsystems A and B, respectively. One can use the clarified structure 
of the density matrix set as the union of orbits obtained by action of the 
unitary group on projectors of rank one with matrix form containing 
only one nonzero matrix element. Then the separable density matrix 
of bipartite composite system can be always written as the sum of 
nin2 tensor products (or corresponding mean tensor product), i.e., in 
(64) the factors are state projectors. Each of tensor products contains 
random unitary matrices of local transforms of the fixed local projector 
for one subsystem and for the second subsystem. It means that an 
arbitrary projector of rank one of a subsystem can be always presented 
in the product form p^^ = u^^p^it^^^, where u"^ is a unitary local 
transform and pA is a fixed projector. 

There are several criteria for the system to be separable. We suggest 
in the next sections a new approach to the problem of separability 
and entanglement based on the tomographic probability description of 
quantum states. The states which cannot be represented in the form 
(63) by definition are called entangled states [38]. Thus the states are 
entangled if in formula (63) at least one coefficient (or more) pi is 
negative which means that the positive ones can take values greater 
than unity. 

Let us discuss the condition for the system state to be separable. 
According to the partial transpose criterion [41], the system is sepa- 
rable if the partial transpose of the matrix pAS (63) gives a positive 
density matrix. This condition is necessary but not sufficient. Let us 
discuss this condition within the framework of positive-map matrix 
representation. For example, for a spin-1/2 bipartite system, we have 
shown that the map of a density matrix onto its transpose belongs 
to the matrix semigroup of matrices C. One should point out that 
this map cannot be obtained by means of averaging with all positive 
probability distributions p]^. On the other hand, it is obvious that the 
generic criterion, which contains the Peres criterion as a partial case, 
can be formulated as follows. 

Let us map the density matrix pAB of a bipartite system onto vector 
PAB- Let the vector pab be acted upon by an arbitrary matrix, which 
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represents the positive maps in subsystems A and B. Thus we get a 
new vector 

Pab = {j1^a®Cb)pab. (65) 

Let us construct the density matrix using the inverse map of the 
vectors onto matrices. If the initial density matrix is separable, the new 
density matrix must be positive (and separable). 

In the case of the bipartite spin- 1/2 system, by choosing Ca = 1 and 
with Cb being the matrix coinciding with the matrix g"^, we obtain 
the Peres criterion as a partial case of the criterion of separability 
formulated above. Thus, our criterion means that the separable matrix 
keeps positivity under the action of the tensor product of two semi- 
groups. In the case of the bipartite spin-1/2 system, the 16x16 matrix 
of the semigroup tensor product of positive contractive maps (52) is 
determined by 24 parameters. Among these parameters, one can have 
some correlations. 

Let us discuss the positive map (52) which is determined by the 
semigroup for the n-dimensional system. It can be realized also as 
follows. 

The nxn Hermitian generic matrix p can be mapped onto essentially 
real n^-vcctor p by the map described above. The complex vector p is 
mapped onto the real vector p^- by multiplying by the unitary matrix 
S, i.e., 

Pr = Sp, p = S-^p,. (66) 
The matrix S is composed from n unity blocks and the blocks 

where j corresponds to a column and k corresponds to a row in the 
matrix p. 

For example, in the case n = 2, one has the vector pj- of the form 



One has the equalities 



A/2Re/9i2 
\/2Im pi2 
\ P22 ) 



f,=f = l^p\ (69) 



(68) 



The semigroup preserves the trace of the density matrix. Also the 
discrete transforms, which are described by the matrix with diagonal 
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matrix blocks of the form 



V 



/ 1 0\ 
10 

0-10 

Vo 1/ 



(70) 



preserve positivity of the density matrix. 

For the spin case, the semigroup contains 12 parameters. 

Thus, the direct product of the semigroup (52) and the discrete 
group of the transform D defines positive map preserving positivity 
of the density operator. One can include also all the matrices which 
correspond to other not completely positive maps. The considered rep- 
resentation contains only real vectors and their real positive transforms. 
This means that one can construct representation of semigroup of 
positive maps by real matrices. 



5. Symbols, star-product and entanglement 

In this section, we describe how entangled states and separable states 
can be studied using properties of symbols and density operators of 
different kinds, e.g., from the viewpoint of the Wigner function or 
tomogram. The general scheme of constructing the operator symbols is 
as follows [36]. 

Given a Hilbert space H and an operator A acting on this space, 
let us suppose that we have a set of operators (j (x) acting transitively 

on H parametrized by ?i-dimensional vectors x = (^i, • • • ^ •^n 

). We 

construct the c- number function /^(x) (we call it the symbol of the 
operator A) using the definition 

/l(x) = IV [aU{^)] . (71) 

Let us suppose that relation (71) has an inverse, i.e., there exists a set 
of operators -D(x) acting on the Hilbert space such that 

A = j /^(x)L>(x) c?x, = / /i(x)1V^(x) dx. (72) 

One needs a measure in x to define the integral in above formulae. 
Then, we will consider relations (71) and (72) as relations determin- 
ing the invertible map from the operator A onto the function /^(x). 
Multiplying both sides of Eq. (2) by the operator tj (x') and taking the 
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trace, one can satisfy the consistency condition for the operators U{x.') 
and -D(x) 

(73) 



Tr 



?7(x')i?(x)J = (5 (x' - x) . 
The consistency condition (73) follows from the relation 

/^(x) = J i^(x,x')/^(x')cix'. 



(74) 



The kernel in (74) is equal to the standard Dirac delta-function, if the 
set of functions /^(x) is a complete set. 

In fact, we could consider relations of the form 



and 



^--/i(x) 



/i(x) - A. 



(75) 



(76) 



The most important property of the map is the existence of the asso- 
ciative product (star-product) of the functions. 

We introduce the product (star-product) of two functions /^(x) and 
(x) corresponding to two operators A and B by the relationships 



/iB W = /i(x) * /b(x) := Tr i^J7(x) 



(77) 



Since the standard product of operators on a Hilbert space is an asso- 
ciative product, i.e., A{BC) = {AB)C, it is obvious that formula (77) 
defines an associative product for the functions /^(x), i.e., 

/^(x) * (/^(x) * /^(x)) = (/^(x) * /^(x)) * /^(x). (78) 

Using formulae (71) and (72), one can write down a composition rule 
for two symbols /^(x) and /^(x), which determines the star-product 
of these symbols. The composition rule is described by the formula 



/^(x)*/^(x) = J /^(x")/s(x')i^(x",x',x)dx'a!x". 



(79) 



The kernel in the integral of (79) is determined by the trace of the 
product of the basic operators, which we use to construct the map 



K(x",x',x) = Tr £>(x")-D(x')C>(x) 



(80) 



The kernel function satisfies the composition property K * K = K . 



MMSZ-Bregenz.tex; 1/02/2008; 8:18; p. 23 



24 



V.I. Man'ko, G. Majrmo, E.C.G. Sudcirshan and F. Zaccaria 

6. Tomographic representation 



In this section, we will consider an example of the probability repre- 
sentation of quantum mechanics [42] . In the probability representation 
of quantum mechanics, the state is described by a family of probabil- 
ities [43—45]. According to the general scheme, one can introduce for 
the operator A the function /^(x), where 

X = {xi,X2,X3) = I/), 

which we denote here as 'w^{X, /j,,!^) depending on the position X and 
the parameters /x and v of the reference frame 

w^{X, /i, z/) = TV Uc/(x)] . (81) 

We call the function w^{X, the tomographic symbol of the oper- 
ator A. The operator U (x) is given by 

[/(x) = UiX,ii,u) = exp (^y (qp + pq)^ exp (^^ (g^ + p^)^ X){X \ 

X exp {-^-^ [(f +P^)^ exp ^-y {qp + pq)^ 

= U^,\X){X\UI. (82) 

The tomographic symbol is the homogeneous version of the Moyal 
phase-space density. The angle 9 and parameter A in terms of the 
reference phase-space frame parameters are given by 

jU = e'^cos^, z/ = e~^sin^, 

that is, q and p are position and momentum operators 

q\X)=X\X) (83) 

and I I is the projection density. One has the canonical transform 

of quadratures 

X = U^vq = IJ'q + i^P, 

I^-U^,pU^,- 2^ p 2^ q. 

Using the approach of [46] one obtains the relationship 
[/ {X, /X, v) = 6{X — fiq — up). 
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In the case we are considering, the inverse transform determining the 
operator in terms of the tomogram [see Eq. (72)] will be of the form 

A = j wj^{X,iJL,v)b{X,iJi,v)dXdiJidv, (84) 

where ^ 

£)(x) = D{X, H^f) = — exp {iX — ivp — ijiq) . (85) 
27r 

The trace of the above operator, which provides the kernel deter- 
mining the trace of an arbitrary operator in the tomographic represen- 
tation, reads 

TV£>(x) = e^^(5(/i)(5(z/). 
The function 'u;^(X, /x, z^) satisfies the relation 

wj^ (AX, A/i, \v) = 1^ w^{X, /X, u). (86) 

This means that the tomographic symbols of operators are homoge- 
neous functions of three variables. 

If one takes two operators Ai and A2, which are expressed through 
the corresponding functions by the formulas 

Ai = J w^^{X',fi',i'')b{X',n',u')dX'dn'di^', 

(87) 

A2 = J w^^{X",^",u")b(X",^",u")dX"dfi"diy", 

and A denotes the product of Ai and A2, then the function w^{X, n, u), 
which corresponds to A, is the star-product of the functions w^^ (X, /x, v) 
and w^^{X, Hjv). Thus this product 

{X, n,v)= w^^ {X, w^^ {X, n, u) 

reads 

w^{X,fi,ij) = J w;^^(x")u;^^(x)i^(x",x',x)dx"dx', (88) 
with kernel given by 

K{x", x', x) = Tr [b{X", n", v")b{X' , u')U{X, /x, i/)l . (89) 

The explicit form of the kernel reads 

K{Xi,lli,Vi,X2,lJL2,^2,X^,v) 
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6 (//(i/i + 1/2) - i^im + 1^2) 



47r2 



exp f ^1 {uiH2 - + 2X1 + 2X2 



[1^1 + 1^2) H (a^i + fJ-2) 



X 



(90) 



7. Multipartite systems 

Let us assume that for multimode (A^-mode) system one has 



U{y)=l[m{s^''^), (91) 

k=l 

N 

b{y) = Jl(E)D (f^'^)) , (92) 



k=l 



where 



This means that the symbol of the density operator of the composite 
system reads 

N 

fp{y) = T:v[p^®U{x^^))\. (94) 
fe=i 

The inverse transform reads 

N N m 



/J Y 1 y I II, 

dyfM n ®Dix(% = n n ^4'^- 
k=l k=ls=l 



(95) 



Now we formulate the properties of the symbols in the case of 
entangled and separable states, respectively. 

Given a composite m-partite system with density operator p. 

If the nonnegative operator can be presented in the form of a "prob- 
abilistic sum" 

p = Y.nz)p^^'^ C5 p^P^ C5 • • • C5 4^""), (96) 



with positive probability distribution function V{z)^ where the com- 
ponents of z can be cither discrete or continuous, we call the state a 
"separable state." Without loss of generality, all factors in the tensor 
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products can be considered as projectors of rank one. This means that 
the symbol of the state can be presented in the form 

m 

fpiy) = T.nz)Y[&''\^k,z). (97) 

2 k=l 

Analogous formula can be written for the tomogram of separable state. 
We point out that in the multipartite case one can introduce ran- 
dom symbols and represent the symbol of separable density matrix 
of composite system as mean value of pointwisc products of symbols of 
subsystem density operators. As in the bipartite case, one can use sum 
over different indices but this sum can be always reduced to the sum 
over only one index common for all the subsystems. It is important that 
for separable state its symbol always can be represented as the sum 
containing number of summants which is equal to dimensionality of 
composite system. Each term in the sum is equal to mean value of ran- 
dom projector. The random projector is constructed as the product of 
diagonal projectors of rank one in each subsystem considered in random 
local basis obtained by means of random unitary local transforms. 



8. Spin tomography 



Below we concentrate on bipartite spin systems. 

The tomographic probability (spin tomogram) completely deter- 
mines the density matrix of a spin state. It has been introduced in 
[31, 32, 36]. The tomographic probability for the spin-j state is defined 
via the density matrix by the formula 

{jm I D\g)pD{g) \ jm) = W^^\m, 0), m = -j, -j + 1, . . . , j, 

(98) 

where D{g) is the matrix of S'?7(2)-group representation depending on 
the group element g determined by three Euler angles. It is useful to 
generalize the construction of spin tomogram. 

One can introduce unitary spin tomograms w{m, u) by replacing in 
above formula (98) the matrix D{g) by generic unitary matrix u. For 
the case of higher spins j = 1, 3/2, 2, . . ., the nxn projector matrix 



Pi 



/I 


• 


■ °\ 





• 


• 


Vo 


• 


• 0/ 



n 



2j + l 



(99) 
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has the unitary spin tomogram denoted as 

Wl{j,U) = \uu\'^, Wi{j - l,u) = {Uul"^, ... Wi{-j,U) = {Uinl"^. 

(100) 

Other projectors 

/O ••• 



Pk 







\0 • •• 





oj 



(101) 



in which unity is located in kth column, have the tomogram W};(m,u) 
of the form 

'Wk{j,u) = \Uklf, Wk{j -1,U) = \Uk2f, ... Wk{-j,u) = \Ukn\'^. 

(102) 

In connection with the decomposition of any density matrix in the form 



P — ^Pjk^jk, 
jk 



(103) 



the unitary spin tomogram can be presented in form of the decompo- 
sition 

Wp{m, u) = Y^ pjkWjk{m, u), (104) 

jk 

where Wjk{rn, u) are basic unitary spin symbols of transition operators 
Ejk of the form 



Wjk{m,u) = {jm I u'EjkU \ jm). 

If one uses a map 

P^ p', 

the unitary spin tomogram is transformed as 

Wp{m,u) w'p{m,u) = ^ Pjf,Wjk{m,u). 

jk 

If the transform (106) is a linear one 

Pjk ^ Pjk ~ ^jk,psPpsi 

the transform reads 

w'p{m, u) = Y^ ppsw'pg{m, u). 

ps 



(105) 
(106) 

(107) 

(108) 
(109) 
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Here 

w'psim, u) = Y^ Ljk,psWjk{m, u) (110) 

jk 

is the hnear transform of the basic tomographic symbols of the opera- 
tors Ejj^. 

Let us now discuss some properties of usual spin tomograms. 

The sot of the tomogram values for each is an over complete set. 
We need only a finite number of independent locations which will give 
information on the density matrix of the spin state. Due to the structure 
of the formula, there are only two Euler angles involved. They are 
combined into the unit vector 

= (cos (?!) sin 7?, sin ^ sin??, COST?). (HI) 

This is the Hopf map from to 5^. 

The physical meaning of the probability W{m,0) is the following. 

It is the probability to find, in the state with the density matrix p, 
the spin projection on direction equal to m. For a bipartite system, 
the spin tomogram is defined as follows: 

W{mim20i02) = (ji"iii2"i2 I D\gi)D\g2)pD{gi)D{g2) \ jimij2m2). 

(112) 

It completely determines the density matrix p. It has the meaning of the 
joint probability distribution for spin ji and j2 projections mi and m2 
on directions Oi and O2. Since the map p ^ W is linear and invertible, 
the definition of separable system can be rewritten in the following form 
for the decomposition of the joint probability into a sum of products 
(of factorized probabilities): 

W{mim26i62) = ^PkW'^''\mi6i)W^''\m202)- (113) 

k 

This form can be considered to formulate the criterion of separability 
of the two-spin state. 

One can present this formula in the form 

^(mimaOiOs) = (W"(mi0i)#(m202)), (114) 

where we interpret the positive numbers as probability distributions. 
Thus separability means the possibility to represent joint probability 
distribution in the form of average product of two random probability 
distributions. 

The state is separable iff the tomogram can be written in the form 
(113) with X^fcPfc = 1) Pfc ^ 0. It seems that we simply use the definition 
but, in fact, we cast the problem of separability into the form of the 
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property of the positive joint probability distribution of two random 
variables. This is an area of probability theory and one can use the 
results and theorems on joint probability distributions. If one does not 
use any theorem, one has to study the solvability of relation (113) 
considered as the equation for unknown probability distribution and 
unknown probability functions ^^('^^(miOi) and T^('^)(m202). 



9. Example of spin-1/2 bipartite system 

For the spin-1/2 state, the generic density matrix can be presented in 
the form ^ 

p=-{l + a-n), n = (ni, 7x2,7x3), (115) 

where a are Pauli matrices and fp < 1, with the vector n for a pure 
state being the unit vector. This decomposition means that we use as 
basis in 4-dimensional vector space the vectors corresponding to the 
Pauli matrices and the unit matrix, i.e.. 



0-1 



1 
1 

\0j 



0-2 



/ \ 

— I 



\ / 



The density matrix vector 



0-3 



/ 1 \ 




V-1/ 





Vi/ 



(116) 



fpii\ 
P12 
P21 

V P22 J 

is decomposed in terms of the basis vectors 

p=^(l+ + 7Z2<T2 + ^3(73) . 



(117) 



(118) 



This means that the spin tomogram of the spin-1/2 state can be given 
in the form 



1 n-0 

2 + ^ 



n- 



(119) 



We can consider tomograms of specific spin state. If the state is pure 
state with density matrix 



P+ 



1 




(120) 
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the spin tomogram W{m,0), where 

m = = (sin ^ cos sin ^ sin cos 0) 



has the values 



W+( J,o) =cos2?, o) =sin2^. (121) 



The tomogram of the pure state 



= ( 2 ; ) , (122) 



has the values 



W-{-,0] = sin^- = cos 



2' 7 2 2 



^_ - = cos^ - = sin^ 



(123) 



2' y 2 2 

The spin tomogram of the diagonal density matrix 
equals 

Wd{m,0) = pnW+{m,0+) + p22W-{mX)- (125) 

The generic density matrix which has eigenvalues pn and P22 can be 
presented in the form 

uoPdul, (126) 

where the unitary matrix uq has columns containing components of 
normalized eigenvectors of the density matrix p. 

This means that the tomogram of the state with the matrix p reads 

Wp{m,0) = {m I u^uopdUQU \ m). (127) 

The elements of the group can be combined 

u^uq = u. (128) 

Thus the tomogram becomes 

Wp{m,0) = Wd{m,0'), (129) 
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where the angle 0' corresponds to the Euler angle calculated from the 
product of two unitary matrices UqU. 

One can use the property of numbers pu and P22 to interpret formula 
(125) as averaging 

Wd{m,6) = {Wim,6')), (130) 

where one interprets two functions TV+(m, 0) and W-{m,u) as the 
realization of "random" probability distribution function W±{m,0). 
Then one has 

Wp{m, 0) = {W{m, 0')). (131) 

The new vector 0' has the parameter 9' obtained from the initial pa- 
rameter 9 by action of the unitary matrix on the initial unitary matrix 
u. 

Inserting these probability values into relation (113) for each value 
of fc, we get the relationships 




k 



w 01,02) = \ + l (Ef^^^fe) • oi - ^ i^P'^*^ ■ 

-J2Pki^k-0i){nl-62), (133) 



(-^,^,01,02) = l-i|^^p,nfcj -Oi + i^^X^p.n^j -02 

-J2Pk{^k-0i)(ril-62). (134) 



k 

One has the normalization property 



1/2 

J2 W^(TOim20i02) = 1. (135) 

mi,m2=— 1/2 



One easily gets 



W^Q,^,0i,02) +T^Q,-^, 01,02) + (^^Pfcnfcj -Oi. (136) 
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This means that the derivative in Oi on the left-hand side gives 



W 



d 

Analogously 

d 

9^ 



1 1 
2' 2 



,0i,02 




w 



+ W {—,^,61,02 




Taking the sum of (133) and (134)) one sees that 
1 d d 



2 dOi dOj 



+ w 



2'2'°^'°^ 



n 



(137) 



(138) 



(139) 



Since we look for the solution where pk > 0, we can introduce 



Pknk, 



iV, 



-,(*) 

PkH'. 



(140) 



This means that the derivative in (139) can be presented as a tensor 

-T,j = Y.iNkUNi%. (141) 



(142) 
(143) 



One has 

'^Pkfik = J2VPk^k, 
k k 

j:pk^k=j:v^kNk^- 

k k 

The conditions of solvability of the obtained equations is a criterion for 
separability or entanglement of a bipartite quantum spin state. Using 
the arguments on the representation of the tomogram (tomographic 
symbol) as sum of random basic projector symbols we get that for two 
qubits the separable state has the tomogram with following properties. 
All four values of joint probability distribution function are equal to 
mean values of product of two cosine of two different angles squared, 
product of sine of two different angles squared and product of sine and 
cosine squared, respectively. The entangled matrix does not provide 
such structure. 

As an example, we consider the Werner state. For the Werner state 
(see, e.g., [47]) with the density matrix 

2 \ 



0^0 



I 1+2 
' 4 



PAB 



\ 



4 u 



^-^j 
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(144) 



PA 



PB = ^ 



1 
1 



one can reconstruct the known results that for p < 1/3 the state is sepa- 
rable and for p > 1/3 the state is entangled, since in the decomposition 
of the density operator in the form (113) the state 



Po 



/ 1 0\ 

10 

10 

\0 1/ 



(145) 



has the weight pq = {1 — 3p)/4. 

For p> 1/3, the coefficient Po becomes negative. 

There is some extension of the presented consideration. 

Let us consider the state with the density matrix (nonnegative and 
Hermitian) 



Ri2\ 

Pll Pl2 

P21 P22 

Vi?21 R22J 



Tip= 1. 



(146) 



Using the procedure of mapping the matrix onto vector p and applying 
to the vector the nonlocal linear transform corresponding to the Peres 
partial transpose and making the inverse map of the transformed vector 
onto the matrix, we obtain 



/ pu \ 

Pll R12 

R21 P22 

Vp21 1^22/ 



(147) 



In the case of separable matrix p, the matrix p"^ is a nonnegative 
matrix. Calculating the eigenvalues of p^ and applying the condition 
of their positivity, we get 



-R1I-R22 > |P12P, PIIP22 > I-R12I' 



(148) 



Violation of these inequalities gives a signal that p is entangled. For 
Werner state (144), Eq. (148) means 

l+p>0, l-p>2p, (149) 

which recovers the condition of separability p < 1/3 mentioned above. 
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The joint probability distribution (112) of separable state is positive 
after making the local and nonlocal (partial transpose-like) transforms 
connected with positive map semigroup. But for entangled state, func- 
tion (112) can take negative values after making this map in the func- 
tion and replacing on the right-hand side of this equality the product of 
two matrices D{g) by generic unitary transform u. This is a criterion 
of entanglement in terms of unitary spin tomogram of the state of 
multiparticle system. 

A simpler and more transparent case is the generalized Werner 
model with density matrix 

P=^{1+ lliai (g) Ti -I- /X20-2 (8) T2 + H3(^3 ® T3) • (150) 

Here the density matrix is expressed in terms of tensor products of two 
sets of Pauli matrices and {k = 1,2,3), which are chosen in the 
standard form. 

Its eigenvalues are 

l-Hl-Hl-HZ, l + Atl+A*2-/^3, l + Atl-/^2 + /^3, I-/X1+/X2 + /X3. 

These eigenvalues are related to the vertices of a regular tetrahedron. 
The partially time-reversed density matrix is 

^ = ^ (1 - /ilO"! (X" n - /X2CJ2 ®T2- 1^3(^3 <E) T3) , (151) 

which may be viewed as 

P = l('^^l(''^P with LW/,W=pW, l(V^ = 1-/>(2). 

(152) 

The eigenvalues of this are 

I + //I+//2 + /X3, l + /Ui-/V,2-/i3, l-/ii+/i2-/i3, I- Hi- 112 + ^3- 

These form an inverted tetrahedron and they have the common domain 
which is a regular octahedron. The unitary spin tomograms can be 
written down by inspection and we may verify that all the relations 
required by the separability criterion (see the next section for details) 
are satisfied by any point inside the octahedron for p and for L^^^^L^^^p 
but the relations connected with positivity condition expressed in terms 
of positivity of unitary spin tomogram fail when it lies outside. 
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10. Tomogram of the group U{n) 

In this section we discuss in more detail the separabihty criterion using 
introduced notion of unitary spin tomogram. 

In order to formulate a criterion of separability for a bipartite spin 
system with spin ji and j2, we introduce the tomogram w(l,m, g^"'^) 
for the group U (n) , where 

n = nin2, ni = 2ji + 1, n2 = 2j2 + 1, 

and g^"^") are parameters of the group clement. Vectors I and m label a 
basis I Lm) of the fundamental representation of the group U{n). For 
example, since this representation is irreducible, being reduced to the 
representation of the U{ni) (g) U{n2) subgroup of the group U{n), the 
basis can be chosen as the product of basis vectors: 

\ji,mi) \ j2,m2) =\ ji,j2,mi,m2). (153) 

Due to the irreducibility of this representation of the group U (n) and 

its subgroup, there exists a unitary transform u^j^j^rnim2 I ^''^) 
that 

I ii,j2,mi,m2) = Xl^^rj2mim2 I ^"^). (1^4) 

lm 

I M = (""^)S2mim2 I jlJ2,mi,m2). (155) 

mim2 

One can define the U (n) tomogram for a Hermitian nonnegative nxn 
density matrix p, which belongs to the Lie algebra of the group U (n) , 
by a generic formula 

w{l,m,g^"'^) = {T,in\ U\g^"'^)pU{g^"-^) \l,m). (156) 

Formula (156) defines the tomogram in the basis | I, m) for arbitrary 
irreducible representation of the unitary group. But below we focus 
only on tomograms connected with spins. 

Let us define the U{n) tomogram using the basis | ji, j2,mi,m2) 
namely for fundamental representation, i.e., 

= (Ji,j2,mi,m2 I ?7^(ff("))pf/(5^")) | ii, i2, mi, ms). (157) 

This unitary spin tomogram becomes the spin-tomogram [34] for the 
g{n) G [/(2) (g) U{2) subgroup of the group U{n). The properties of this 
tomogram follow from its definition as the joint probability distribution 
of two random spin projections mi,m2 depending on ^f^**) parameters. 
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One has the normahzation condition 

^ w;(^i'^2)(mi,m2,5(")) = 1. (158) 

mi,m2 

Also all the probabilities are nonnegative, i.e., 

u;(^'i'J2)(mi,m2,5("^) > 0. (159) 

Due to this, one has 

J2 \w^''^^Hmi,m2,g^''^)\ = l. (160) 

mi ,m,2 

For the spin-tomogram, 

3(^)^(01,02) (161) 

and 

u;(^i'^2)(mi,m2,5(")) ^ w{mi,m2, 01,62). (162) 

The separability and entanglement condition discussed in the previ- 
ous section for a bipartite spin-tomogram can be considered also from 
the viewpoint of the properties of a C/ (n) tomogram. If the two-spin 
nxn density matrix p is separable, it remains separable under the action 
of the generic positive map of the subsystem density matrices. This map 
can be described as follows. 

Let p be mapped onto vector p with components. The components 
are simply ordered rows of the matrix p, i.e., 

P= (P11,P12, • • • ,Pln,P21,P22, • • • ,Pnn,)- (163) 

Let the n^xn^ matrix L be taken in the form 

L = J2ps4''^^4''\ Ps>0, Y1p^ = 1, (164) 

s s 

where the nixni matrix Li^^-* and the n2Xn2 matrix lI"'^^ describe the 
positive maps of density matrices of spin-ji and spin-j2 subsystems, 
respectively. We map vector p onto vector pL 

PL = Lp (165) 

and construct the nxn matrix p^, which corresponds to the vector pL. 
Then we consider the U{n) tomogram of the matrix pi, i.e., 

w^i''''\mi,m2,g^^^) 

= {ji,j2,mi,m2 I ?7t(5("))pi?7(5W) | j^,j2,mi,m2). (166) 
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Using this tomogram we introduce the function 

F{g^^\L)= |4^'^'^'^)(mi,m2,5("))|. (167) 

mi,m2 

For separable states, this function does not depend on the U (n)-group 
parameter g^"^^ and positive-map matrix elements of the matrix L. 

For the normalized density matrix p of the bipartite spin-system, 
this function reads 

F{g^^\L) = l. (168) 

For entangled states, this function depends on g^"^^ and L and is not 
equal to unity. This property can be chosen as a necessary and sufficient 
condition for separability of bipartite spin-states. We introduce also 
tomographic purity parameter /x^ of feth order by the formula 

For identity semigroup element L and specific g^''^ unitary transform 
diagonalizing the density matrix, the tomographic purity /X2 is identical 
to purity parameter of the state p. The parameters for k = 2, 3, . . ., 
correspond to Trp'^"^^. 

In fact, the formulated approach can be extended to multipartite 
systems too. The generalization is as follows. 

Given N spin-systems with spins ji,j2, ■ ■ ■ ,jN, let us consider the 
group U{n) with 

N 

'^=11'^^' '^k = "^jk + 1- (169) 

k=l 

Let us introduce the basis 

N 

k=l 

in the linear space of the fundamental representation of the group U{n). 
We define now the U{n) tomogram of a state with the nxn matrix p: 

Wp{m,g^''^) = {m \ U\g^''^)pU{g^"^) \ m). (171) 

For a positive Hermitian matrix p with Tr p = 1, we formulate the 

criterion of separability as follows. 

Let the map matrix L be of the form 

N 

^ = E^'^(n^4'0' ps>^^ 12ps = ^^ (172) 

s k=l s 
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where Ls is the positive-map matrix of the density matrix of the 
A;th spin subsystem. We construct the matrix as in the case of the 
bipartite system using the matrix L. The function 

F{g(^\L)=J2H,{m,g^^^)\>l (173) 

m 

is equal to unity for separable state and depends on the matrix L and 
U (n)-parameters g^"^^ for entangled states. 

This criterion can be applied also in the case of continuous variables, 
e.g., for Gaussian states of photons. Function (173) can provide the 
measure of entanglement. Thus one can use the maximum value (or a 
mean value) of this function as a characteristic of entanglement. In the 
previous section, we considered the generalized Werner states. Using 
the above criterion, one can get the domain of values of the parameters 
of the states for which one has separability or entanglement. In fact, 
the separability criterion is related to the following positivity criterion 
of finite or infinite (trace class) matrix A. The matrix A is positive iff 
the sum of moduli of diagonal matrix elements of the matrix UAW 
is equal to a positive trace of the matrix A for an arbitrary unitary 
matrix U. 



11. Dynamical map and purification 

In this section, we consider the connection of positive maps with pu- 
rification procedure. In fact, formula 

p^p' = J2PkUkpUl (174) 

k 

where Uk are unitary operators, can be considered in the form 

P^p' = ^PkPk, Pk>0, ^pk = l. (175) 

k k 

Here the density operators p^ read 

Pk = UkpUl (176) 

and the maps which are not sufficiently general keep the most degener- 
ate density matrix fixed. This form is the form of probabilistic addition. 
This mixture of density operators can be purified with the help of a 
fiducial rank one projector Pq 

p' -^p" = N 



Y^VPkPj 



kj 



PkPpPj 

y/Ti pkPoPjPo 



(177) 
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where iV is a normalization constant 

The normahzation is unnecessary if all pk arc mutually orthogonal. We 
call this map a purification map. It maps the density matrix of mixed 
state on the density matrix of pure state. 

The map (174) could be interpreted as the evolution in time of the 
initial matrix po considering unitary operators Uk{t) depending on time. 
Thus one has 

Po ^ p{t) = Y.PkUk{t)poUl{t). (179) 

k 

In this case, the purification procedure provides the dynamical map of 
a pure state 

I V'o)(^o HI V'(i))(V'(i) I, (180) 

where | ip{t)) obeys a nonlinear equation and, in the general case, this 
map does not define a one parameter group of transformations not even 
locally. 

For some specific cases, the evolution (179) can be described by a 
semigroup. The density matrix (179) obeys then a first-order differential 
equation in time for this case [27-29] . 

More specifically, the reason why there is no differential equation in 
time for the generic case is due to the absence of the property 

Pu(i2) = Y.K^''{t2,h)Pmn{tl), (181) 

mn 

where the kernel of evolution operator satisfies 

K^"{t^,t2)K^^{h,h) = Kff{ts,h). (182) 

Thus, via a purification procedure and a dynamical map applied to 
a density matrix we get a pure state (nonlinear dynamical map). This 
map can be Tiscd in nonlinear models of quantum evolution. Many linear 
positive maps both completely positive and not completely positive are 
contractive. We define a positive map L as "contractive" or "dilating" 
if Tr {Lp)"^ < Tr (p)^ or Tr (Lp)'^ > Tr (p)^, respectively. This means, for 
example, that purity parameter p = Ti p^ after performing the positive 
map generically becomes smaller. There are maps for which the purity 
parameter is preserved, for example, 

p^ptr^ p^-p+^1. (183) 



MMSZ-Bregenz.tex; 1/02/2008; 8:18; p. 40 



The geometry of density states, positive maps and tomograms 41 

These linear maps include also unitary transform 

p-^upv). (184) 
There are maps which provide dilation. For qubit system, matrix 



L = 



/ 1 1 \ 



Vo 0/ 



(185) 



acting on arbitrary vector pb corresponding to a density matrix po gives 
the pure state 

"'-{ID- 

The matrices the inverse matrices exist for e 7^ 0. But these inverse 
matrices do not provide positive trace preserving maps. Since the purifi- 
cation procedure provides a positive map, which increases the purity 
parameter, the composition of linear map with the purification map 
provides the possibility to recover the initial density matrix p which 
was the object of action of a positive linear map. It means that the 
purification map Lp can give 

Lp,{Lp) = lp (187) 

for any density matrix p but the choice of fiducial projector depends 
on p (the initial condition). 

Thus one has also for completely positive maps 



'-T.P'k^ Pk = VkpVl ^F,t^, = l. (188) 



P 

k 



Making polar decomposition 

Pk = y/POkUk, UkUl = 1, pok > 

and introducing the positive numbers pk = Trpofc, we construct the 
map 



P' 



n 1 , PkPoPj + PjPoPk I 1 

P =<2^VPkPj /Tn PnP r ' 2^Pk = ^. PkPk = Pk- 

[ kj VTpknpjPo J 

(189) 

The matrix p" is a matrix of rank one for any rank one fiducial projector 

Pq. The projector Pq is restricted to be not orthogonal to the generic 
matrix pk- Taking A'' orthogonal projectors Pq^-* (s = 1, 2, . . . , A'') and 
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obtaining N projectors p", one can combine them in order to get the 
initial matrix p. It means that one can take convex sum of the N 
pure states p" to recover the initial mixed state p. Another way to 
make the state with higher purity was demonstrated using the modified 
purification procedure in [48]. For qubit state, one has 

. P1P0P2 + P2P0P1 . . 

P = PiPi + P2P2 + kVPiP2 ^y^YFpiPop2Po ' ^^1+^^2 = 1, (190) 

where the decoherence parameter < «; < 1 is used. If k ~ 1, we 
increase purity. 

Let us discuss the map (188) using its matrix form, i.e., 

Pa/3 ^ Pap = Yl ^c.fS,ijPij- (191) 

The matrix J0.a/3,ij is expressed in terms of the matrices Vk as 

= J2(^k)a^mpj. (192) 
k 

One can construct another positive map [49] 

p^p' = ^rfeTV(i?fep), (193) 

k 

where are density matrices and are positive operators satisfying 
the normalization condition 

k 

The matrix corresponding to this map (called entanglement breaking 
map [50]) reads 

^ip,ij=J2i^k)a/3iR*k)iJ. (195) 

k 

The entanglement breaking map is contractive positive map. There 
exist some special cases of completely positive maps. For example, 

P^-ep + ^p (196) 

differs from the depolarizing map by replacing the unity operator by 
the density operator. Another map reads 
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The decoherence map (phase damping map) of the kind 

Pij^l f'' 'Z^- (198) 

where |A| < 1 provides contractive map with uniform change of off- 
diagonal matrix elements of the density matrix. 

Let us discuss the property of tomogram of bipartite system with 
density matrix pi2- If the density matrix is separable, than the depo- 
larizing map of the second subsystem provides the following density 
matrix 

P12 ^Pe = -£P12 + P^^h, (199) 

where 

p«=Tr2(pi2) (200) 

and I2 is the A''2-dimensional unity matrix. Then one has the property 
of unitary spin tomogram 

We{mi,m2,g^''^) = -£Wi2{rni,m2,g^''^)+^^^w{mi, 1712,9^"-^), (201) 

iV2 

where ^f^"-* is matrix of u(^{2ji + l)(2j2 -|- 1)^ unitary transform; 

We{mi, 1712,9^"^^) is the tomogram of transformed density matrix of bi- 
partite system; 

w{Tni,m2, g^"^^) is the unitary spin tomogram of tensor product of par- 
tial trace p^^-* over the second subsystem's coordinates of the density 
matrix pi2 and unity operator I2; 

1^12 (wii , m2 , 5^"^ ) is the unitary spin tomogram of the state with density 
matrix pi2. 

The criterion of separability means 



1 + e 



E E 

mi = -jl 1712= -j2 



2j2 + 1 



w{mi,m2,g^'^^) - ewi2 (mi , m2 , 5^"^ ) 



= 1 
(202) 

for arbitrary g^'^^ and e. 

For Werner states pyy, the tomogram of transformed state (in this 
case, it means that p —sp) is related to the initial-state tomogram 
WW 

We (mi , m2 , 5^"^ ) = -£Wi2 (mi , m2 , 5^"^ ) + — ;— • (203) 



The criterion of separability yields 

1. (204) 



1/2 

— eww{mi,m2,g'^'"'') 



E 

mi,m2=— 1/2 
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Equality (204) takes place for arbitrary g^'^^ and e only for \p\ < 1/3. 
For p > 1/3, the above sum depends on g^'^^ and e and it is larger than 
one. 

It is obvious if one calculates the tomogram using the element of the 
unitary group of the form 



in) 
% 



/O 1\ 
10 
10 

V 1 0/ 



(205) 



At this point, the sum (204) reads 



1/2 

E 

mi,m2=- 



1/2 



1+e 



eww{mi,m2,g 



(n)^ 



= 3 


l + ep 


+ 


1 — 3pe 




4 


4 



(206) 

One can see that this sum equals to one independently on the value 
of parameter |e| < 1 only for values \p\ < 1/3. For p = 1, the maxi- 
mum value of the sum equals 2 = (1 + 3£)/2 (e = 1). This value can 
characterize the degree of entanglement of Werner state. 

We have introduced positive nonlinear map of density matrix which 
is purification map. The purification map can be combined with con- 
tractive maps discussed. The tomographic-probability distributions un- 
der discussion can be completely described by their characteristic func- 
tions. This means that the relation of tomogram property to entan- 
glement can be formulated in terms of the properties of characteristic 
functions. 

One can also check the criterion using example of two-qutrite pure 
entangled state with wave function 



V') 



1 

71 



(207) 



m=— 1 



The sum defining the criterion of separability for specific U (9) trans- 
form (/q"^ which is diagonalizing the hermitian matrix \ if)) {if) \ 
reads 

For 1/2 > e > 1/8, this sum is larger than one, that means that the 
state is entangled. For e = 1/2, the function has maximum and it is 
equal to 5/3. 

The entanglement of the considered state can be detected using 
partial transposition criterion too. 
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For the case of pure entangled state of two-qutrite system with the 
wave function 

\^) = ^[\ui)\vi)+\uo)\vo)), (209) 

in which the states with spin projections m = — 1 do not participate, the 
partial transpose criterion docs also detect entanglement. Our criterion 
yields for specific U (9) transform , which diagonalizes the hermitian 
matrix | | the following expression for the function F{£,gl^^), 

which reads 

F(.,sr) = 5ii±£l + ii^. (210) 

The function takes maximum value for e = 1/2 that equals to 3/2. 
This value is smaller than 5/3 of the previous case. It corresponds 
to our intuition that the superposition of three product states of two 
qutrite system is more entangled than the superposition of only two 
such product states. 

The criterion can be extended to multipartite spin system. 

We have to apply for n-partite system the transform of the density 
matrix p of the form 

L,^=L«C5Lg)c5...C5 4«), (211) 

(k) 

where the transform Ls^. acts as depolarizing map on the kth subsys- 
tem. If the state is separable 

k k 

each of the terms p^"''' (j = 1, 2, . . . , n) in the tensor product is replaced 
by the term 

pi^'^---.pl^Vi^l,. (213) 
This means that the transformed density matrix reads 

L^P = ^Pk 
k 

The unitary spin tomogram of the transformed density matrix takes 
the form {£ = £i,£2, . . ■ , £n) 

w^{mi,m2,.. .,mn,g^^^) = Y^PkW^';'' {mi,m2, . . .,mn,g^^\£), 

(215) 



(214) 
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where N = n"=i(2is + 1) and element g^^^ is the unitary matrix in 
A?^-dimensional space. The tomogram w^\mi,m2, ■ ■ ■ ,mn,g^^\s) is 
the joint probabihty distribution of spin projections = —js, —js + 
l,...,js, which depends on the unitary transform the state with 

density matrix 

Pk = fl^[-eJ^^ + ^lsy (216) 

For the elements 

n 

= n ^M'^Js + 1), 

where Us{2js + 1) is unitary matrix, the tomogram (215) takes the form 
of sum of the products 

n 

wg{mi,m2,.. .,mn,glr^) = ^Pk H Wkims,Us{2js + l),es), (217) 

k s=l 

with Wk{ms,Us{2js + l),£s) being the unitary spin tomograms of the 

(s) 

sth spin subsystem with transformed density matrix L^^p), . If one uses 
as the matrix Us{2js + 1), the matrix of unitary irreducible represen- 
tation of the SU (2) group, the tomogram Wk depends only on the two 
parameters defining the point on the sphere S'^. 

For a separable state of the multipartite system, one has 

{NY 



^ U,-(mi,m2,...,m„,5W) =1 (218) 



mi,...,m„ 



for all elements g'^^'^ and all parameters £. 

For entangled state, there can be some values of parameters e and 
group elements g^^'^ for which the sum is larger than one. 



12. Conclusions 

We summarize the results of the paper. 

The notion of entangled states (first discussed by Schrodinger [4, 51]) 
has attracted a lot of efforts to find a criterion and quantitative charac- 
teristics of entanglement. A criterion based on partial transpose trans- 
form of subsystem density matrix (complex conjugation of the subsys- 
tem density matrix or its time reverse) provides the necessary and 
sufficient condition of separability of the system of two qubits and 
qubit-qutritt system [52]. The phase-space representation of the quan- 
tum states and time reverse transform (change of the signs of the 
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subsystem momenta) of the Wigner function in the case of Gaussian 
state was apphed to study the separabihty and entanglement of pho- 
ton states in [13]. Recently it was pointed out that the tomographic 
approach of reconstructing the Wigner function of quantum state [43- 
45] can be developed to consider the positive probability distribution 
(tomogram) as an alternative to density matrix (or wave function) 
because the complete set of tomograms contains the complete informa- 
tion on the quantum state [42] . This representation (called probability 
representation) was constructed also for spin states including a bi- 
partite system of two spins. Up to now the problem of entanglement 
was not discussed in the tomographic representation. Some remarks 
on tomograms and entanglement of photon states in the process of 
Raman scattering were done in [53]. The tomographic approach has 
the advantage of dealing with positive probabilities and one deals with 
standard probability distributions which are positive and normalized. 

We studied the properties of separable and entangled state of mul- 
tipartite system using the tomographic probability distributions. The 
positive and completely positive maps of density matrices [39, 54] in- 
duce specific properties of the tomograms. The properties of the posi- 
tive maps were studied in [55]. We formulated necessary and sufficient 
conditions of separability and entanglement of multipartite systems in 
terms of properties of the quantum tomogram. Since the tomograms 
were shown [36] to be related to the star-product quantization proce- 
dure [56], we discuss entanglement and separability properties in terms 
of generic operator symbols. The tomographic symbols of generic spin 
operators were studied in [36] . Then we focused on properties of entan- 
glement and separability of a bipartite system using spin tomograms 
{SU (2)-tomograms) and tomograms of the ?7(A^)-group. 

The idea of the approach suggested can be summarized as following. 

The positive but not completely positive linear maps of a subsystem 
density matrix do preserve the positivity of separable density matrix 
of the composite system. These maps contain also maps which do not 
preserve the positivity of the initial density matrix of an entangled state 
for the composite system. It means that the set of all linear positive 
maps of the subsystem density matrix (this set is semigroup) creates 
from the initial entangled positive density matrix of composite system 
a set of hermitian matrices including the matrices with negative eigen- 
values. To detect the entanglement we use the tomographic symbols 
of the obtained hermitian matrices. The tomographic symbols of state 
density matrices (state tomograms) are standard probabilities. In view 
of this the tomographic symbols of the obtained hermitian matrices 
corresponding to initial separable state preserve all the properties of 
the probability representation including positivity and normalization. 
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But in case of entangled state the tomographic symbols of the obtained 
hermitian matrices can take negative values. The different behaviour 
of tomograms of separable and entangled states of composite systems 
under action of the semigroup of positive maps of the subsystem density 
matrix provides the tomographic criterion of the separability. 

To conclude, we point out the main result of the work. 

We found the criterion of separability which is given by equation (173). 
The criterion is valid for multiparticle spin system. The criterion can 
be called "tomographic criterion" of separability. The tomographic cri- 
terion can be considered also for symplectic tomograms of multimode 
photon states. The condition of separability is sufficient because there 
always exists a unitary group element by means of which any hermitian 
matrix can be diagonalized.lt means that tomographic symbol of non- 
positive hermitian matrix has nonpositive values for some unitary group 
parameters. The suggested criterion is connected with properties of 
the constructed function (173) which for given density matrix depends 
on unitary group parameters g and the parameters of positive map 
semigroup L. For separable density matrix the dependence on unitary 
group parameters and the semigroup parameters disappears and the 
function becomes constant equal to unity. For entangled states the 
function differs from unity and depends on both group and semigroup 
parameters. The suggested criterion can be considered as some com- 
plementary test of separability together with other criteria available in 
the literature (see, for example, [38, 52]). We point out that suggested 
criterion differs from available usual ones by the kind of the necessary 
numerical calculations. To apply this criterion one needs to calculate 
the sum of moduli of diagonal matrix elements of product of three 
matrices. One of the matrices is hermitian and two others are unitary 
ones. This procedure does not need the calculation of the eigenvalues of 
a matrix. The structure of positive (including not completely positive) 
map semigroup with elements L needs extra investigation (see, for 
example, [57]). We found also a test of entanglement based on the 
property of unitary spin tomogram. 

The discussed purification map can be applied to find new quantum 
evolution equations in addition to known ones [58-61]. The application 
of different forms of positive maps [55, 62] and supermatrix represen- 
tation of the maps [63, 64] are useful for better understanding of the 
computations. Entanglement phenomena can be considered using sym- 
bols of density matrix of different kinds, e.g., particular quasidistribu- 
tions [65] as well as tomographic symbols [36] . The difference of symbols 
of entangled and separable density operators for different schemes of 
the star-product quantization needs further investigations as well as 
test of entanglement of some generalizations of Werner state [66, 40] 
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in multipartite case. A relation of tomographic approach to different 
positive maps [67] should be investigated. The tomographic symbols 
are analytic in group parameters. This can be used to find cxtrema of 
tomograms which give information on degree of entanglement. 
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